The experimental observation of quantum Hall effect of I = 3 chiral charge carriers in 

trilayer graphene. 
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Low-energy electronic states in monolayer and bilayer graphenes present chiral charge carriers 
with unique and unusual properties of interest for electronic applications Here, we report the 

magnetotransport measurements in the ABC-stacked trilayer graphene as a function of charge car- 
rier density, magnetic field, and temperature, which show clear evidence of I — 3 chiral quasiparticles 
with cubic dispersion P4lll|. existing in a large, 20 x 60 /im 2 device. Shubnikov-deHaas oscillations 
(SdHO) reveal the Berry's phase 37r, and the marked increase of cyclotron mass near charge neutral- 
ity, consistent with divergent behavior expected for I — 3 quasiparticles. We also observe the pre- 
dicted unconventional sequence of quantum Hall effect (QHE) plateaus, a xy — ±6e 2 //i, ±10e 2 //t, .... 

PACS numbers: 73.43.-f 73.63.-b 71.70.Di 
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The discovery of graphene presented not only realiza- 
tions of truly two-dimensional (2D) electronic systems, 
but revealed entirely new, exotic electronic states. In 
monolayer graphene the low-energy electronic structure 
is described by chiral quasiparticles with linear disper- 
sion, they obey the relativistic Dirac equation and have 
a Berry phase of 7r [2, 0|. In bilayer graphene the 
I = 2 chiral charge carriers have quadratic dispersion 
and Berry's phase 27r. Unusual properties are revealed in 
magnetotransport experiments in a perpendicular mag- 
netic field, where the Berry phase determines the shift 
of Shubnikov-deHaas resistance oscillations, while their 
period is governed by filling of Landau levels (LLs) of 
chiral 2D quasiparticles. The quantum Hall effect shows 
an unconventional sequence of plateaus of Hall conduc- 
tivity, a xy , with quantized steps of 4e 2 //i, except for the 
first plateau, where it is governed by the Berry's phase. 

Despite significant interest in studying layered 
graphene systems with more than two layers, experimen- 
tal progress has been limited (l2T[l6| . Low-energy elec- 
tronic properties depend crucially on the stacking order 
of graphene layers [4-11, [l6j], and therefore such stud- 
ies require samples with well-defined stacking sequence. 
In a bilayer, two honeycomb nets of carbon atoms are 
positioned with half of the atoms of the top layer (B) 
right above the atoms of the bottom layer (A) and the 
other half at the centers of the hexagonal voids in it. The 
third carbon net in a trilayer can either be placed with 
its atoms above the atoms of the bottom layer A, as in 
Bernal structure of crystalline graphite [l7| . or with its 
voids above the lined- up atoms pairs in layers A and B, 
thus breaking the reflection symmetry (Fig. la). The lat- 
ter, ABC stacking, is found in metastable rhombohedral 
modification of graphite (l7| . 

The electronic structure of graphene multi-layers is 
derived from the hybridization of monolayer states via 



interlayer hopping. Its main features are captured al- 
ready by only considering hopping between the nearest- 
neighbour carbons, which are stacked above each other 
in two adjacent layers, 71 ~ O.I70, Fig. la (70 3.16 eV 
is the intra-layer hopping, in bulk graphite 71 0.4 eV, 
and there are also further- neighbour hoppings, 72 — 75, 
which are not shown) [a-lllj. ^ n bilayer, low-energy elec- 
tronic states retain the chiral character but have flatter, 
quadratic dispersion. The effective Hamiltonian is the 
nonlinear generalization of Dirac- Weyl quasiparticles of 
the monolayer, 
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with I = 2. Here, vi = u'/Ti , where v « 10 6 m/s is the 
velocity of Dirac fermions in the monolayer, (p x ,p y ) = 
p(cos</5p,sin<^ p ) is the 2D momentum, tt = p x + ip y , 
a = a x ± iay are the pseudo-spin Pauli matrices, and 
£ = ±1 is a valley index. There are two low-energy val- 
leys with opposite chirality. For our purposes they can be 
considered non-interacting and only yielding a two-fold 
degeneracy of all states; together with the two-fold elec- 
tron spin degeneracy this gives a factor of 4 in the Hall 
conductivity quantization, 4e 2 /h. The wave functions of 
such degree-Z chiral fermions acquire a Berry phase of In 
upon an adiabatic propagation along a closed orbit. In 
bilayer, this results in I = 2, 2tt chiral quasiparticles and 
an unusual integer QHE sequence with a double step, 
a xy = 8e 2 /h, between the hole and electron gases across 
the N=0 LL observed in experiment 

The low-energy band structure of ABA stacked 
graphene trilayer consists of superimposed linear and 
quadratic spectra. Hence, transport is governed by two 
types of chiral quasiparticles: monolayer-like massless 
(/ = 1) and bilayer-like massive (I — 2), albeit with a 
larger effective mass than in bilayer, thaba — V^mAB ~ 
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Figure 1: a, The structure, and b, the low-energy band struc- 
ture H, of ABC-stacked TLG. c, The image of our large-area 
etched Hall bar and schematics of the measurement setup. 
The black scale bar is lOfim. d, Field mobility, n — <j xx / (n a e), 
at T = 0.35 K and e, longitudinal conductivity, as a function 
of BG voltage V g , or carrier density induced via electric field 
effect, n s — C g V g (top scale). The minimum conductivity 
is slightly larger than 6e 2 /h, which agrees with the conduc- 
tance via four quantum channels (dashed line), indicating a 
quantum-mechanical origin. 



0.05m e {niAB = 7i/(2w 2 ) ~ 0.035m e , where m e is the 
electron mass) [5|, HJ. In contrast to bilayer, where appli- 
cation of an electric field across layers opens a band gap 
in the spectrum in ABA trilayer it actually leads to 
a tuneable band overlap fl2- 16|. 

While to our knowledge no magnetotransport experi- 
ments on ABC trilayers have been published, this case 
is actually most interesting, as it is expected to present 
new, I = 3 chiral quasiparticles with cubic dispersion 

GHIllGji, e(p) = 7i0fV7i) 3 > Fig- lb - Tnat this sit- 
uation is quite remarkable could already be seen from 

the fact that the effective mass of such charge carri- 
ers is energy-dependent and diverges at the charge neu- 
trality point (CNP), rriABC = pdp/de = 7 2 /(3pu 3 ) = 
ImAsCTi/e) 1 ^ 3 - This corresponds to a diverging low- 
energy density of states, D{e) ~ e -1 / 3 , in contrast to a 
constant D(e) in bilayer and a vanishing one in mono- 
layer. Such abundance of low-energy scattering states 
would make non-chiral fermions with cubic dispersion 
unstable with respect to decays - but in ABC trilayer 
I = 3 quasiparticles are protected by chirality conserva- 
tion. On the other hand, diverging D(e) leads to super- 
polarizability and super-efficient screening, which were 
predicted theoretically in perpendicular electric field • 
Here we report measurements of magnetotransport os- 
cillations and quantum Hall effect in the ABC-stacked 
trilayer graphene. The optical image of our Hall bar de- 
vice and schematics of the measurement setup are shown 
in Fig. lc. Our sample was prepared following standard 
procedure [20] by mechanical exfoliation of kish graphite. 
While the number of layers can already be identified 
by optical contrast, it was also confirmed using Raman 
micro-spectroscopy, by measuring the full width at half 
maximum (FWHM) of the 2D Raman band [H[23|. The 



Figure 2: Quantum Shubnikov-de Haas Oscillations in ABC 
trilayer. a, Longitudinal resistivity, p X x(V g ), for B = 12 T at 
0.35 K, 4 K, and 10 K, top scale shows carrier density, n s (V g ). 
The inset is a zoom-in of the part surrounded by the dashed 
line, b, SdHO as a function of magnetic field for different 
n a (V g ). The decay of the SdHO magnitude with temperature 
and magnetic field is governed by the cyclotron mass of charge 
carriers, m c (n s ), and by the quantum scattering time, r q . 

characteristic asymmetric shape of the 2D-band spec- 
trum, with dip near 2700 cm" 1 , was used to identify the 
ABC stacking in our sample [23l. l24|. 

The transport properties of the ABC trilayer appear 
enigmatic: a naive estimate of the impurity-limited mo- 
bility following similar arguments to the monolayer case 
gives an energy (charge density) dependent mobility, 
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r 2 / 3 ~ n s . This clearly disagrees with our trans- 
port measurements of the ABC trilayer device of Fig. Ic 
shown in Fig. ld,e. We find that the gate voltage de- 
pendence of longitudinal conductivity, cr xx , is very simi- 
lar to that of monolayer/bilayer graphene, with roughly 
linear asymptotic behavior at high charge carrier densi- 
ties, n s (or gate voltages V g ), for both polarities (+/— 
correspond to electrons/holes). The asymptotic field mo- 
bility estimated assuming a simple Drude model is fi — 
Cxx/(en s ) ~ f 800 cm 2 V _1 s _1 (e is an electron charge) in 
our trilayer, which is lower than in similar high quality 
monolayer samples. Nevertheless, the maximum resis- 
tance near the CNP (V g w 17.3 V), R w 5.7 kfi, roughly 
agrees with the quantum conductance via 4 degenerate 
channels, /i/(4e 2 ) « 6.5 kf2, which is expected for quan- 
tum transport in samples with aspect ratio L/W > f 
(in our case L/W ps f .67). In the opposite limit of wide 
and short samples, one expects minimum conductivity 
(rather than conductance), o™ n ABC — 

The observed monolayer-like "normal" gate voltage de- 
pendence of the conductivity can be reconciled with the 
"abnormal" cubic dispersion of I = 3 chiral fermions if 
one accounts for super-screening of impurity Coulomb 
potential by such quasiparticles. Indeed, in the ran- 
dom phase approximation (RPA), the screened poten- 
tial is governed by the diverging near e = density 
of states, rather than by the bare Coulomb interaction. 
As a result, the divergence in transport scattering cross- 
section is canceled by the vanishing screened potential, 
and one obtains charge-density-independent mobility in 




0.10 0.20 0.25 
f/B (T-<) 



-20 20 40 6' 

V, (V) 
n,(10' z cm- 2 ) 
-6 -4 -2 2 4 



40 




i 


30 




,# 


m 20 




,# 

M 


10 















■ : - - - SG 




: ■ T - - - BG 




: ■■» — TQ 




U o m c 


•v 


V 



20 40 

V,(V) 



Figure 3: a, Landau fan diagram of SdHO. Points for integer 
n correspond to n-th minimum, for half-integers, n — 1/2, to 
n-th maximum of p xx (B) Q, @|- b, Bf(V 9 ) dependence, the 
error bars show standard deviation, the dashed line is linear 
fit. c, The Berry phase of magneto-oscillations from fits in Fig. 
3a. The inset shows quantum scattering time obtained by 
fitting the SdHO by the ALK formula, d, The cyclotron mass, 
m c , of charge carriers as a function of n a (V g ), obtained from 
the same SdHO- ALK fits. Lines show the semi-classical result 
for linear (MG), quadratic (BG) and cubic (TG) spectra, with 
71 = 0.5 eV obtained by fitting the measured m c (n s ). 

the ABC trilayer, ^3 = 8/(7r 2 7Vi). Here Ni is the concen- 
tration of Coulomb impurities, which are mainly located 
in the Si02 substrate. In monolayer graphene screen- 
ing is much weaker, and simply amounts to multiply- 
ing the Coulomb impurity potential by a number, so 
it becomes U(r) = —hva/r. If determined by the di- 
electric constant, £ ~ 4 — 10, of glass substrate alone, 
a = e 2 /(hvT,) ~ 0.5 - 0.2. The mobility of the mono- 
layer is Hi = l/(a 2 TT 2 Ni). Hence, for an equal number of 
impurities, /13 / fi\ — 8a 2 , and the mobility of I = 3 chi- 
ral fermions of the ABC trilayer could appear markedly 
lower than for massless fermions of the monolayer, if a is 
small. 

Having thus established high quality of our ABC 
trilayer sample, we proceeded with magnetotransport 
measurements in a perpendicular field B. Figure 2 
shows examples of SdH quantum magneto-oscillations 
at temperatures 0.35 K, 4 K and 10 K as a func- 
tion of gate voltage for B = 12 T, and as a func- 
tion of magnetic field for several V g . In the semi- 
classical limit of small oscillations, they can be de- 
scribed by the Ando-Lifshitz-Kosevich (ALK) formula 
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sinh x V B 

X = 27r 2 fcsT/ (huj c ), ks is Boltzmann's constant, r q , is 
the quantum scattering time, the cyclotron 

angular velocity, m c = ^r^gj— = ttiabc the cyclotron 
mass, S(e) = irp{e) 2 is the area in the momentum space 
of the orbit at the Fermi energy, Bp = n. s $o/.9LLi where 
$0 = h/e « 4.14-10 -11 T-cm 2 is flux quantum and g^L is 
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Figure 4: QHE and Landau Level spectrum in ABC trilayer 
at T = 0.35 K. a, Longitudinal resistivity, p xx (right scale), 
and Hall conductivity, a xy — p xy /(p xx + p xy ) (left scale) at 
B = 18 T. Vertical arrows with numbers show LL filling, 
v = n s h/ (eB), for the corresponding Quantum Hall states, b, 
o' xy (V g ) for several magnetic fields, B. c, p xx (B) (left scale), 
and p xy (B) (right scale) for V g = 70 V. Quantized plateaus of 
p xy and zero p xx are clearly observed for v — 6, 10. Dashed 
horizontal lines show h/(ve 2 ). Vertical arrows with numbers 
show LL fillings which can be identified in SdHO. d, Schematic 
illustration of the LL DOS in ABC stacked trilayer graphene. 

the LL degeneracy, and j3 — In is the Berry phase of the 
quasiparticles. Fig. 2 already reveals a clear hallmark of 
the large cyclotron mass and small scattering time, r q , 
of I = 3 quasiparticles in the ABC trilayer. The decay 
of the magnitude of SdHO with increasing temperature, 
or decreasing magnetic field, or increasing charge density 
n s , is strikingly fast .much faster than observed in mono- 
layer or bilayer [l|43( (at least half a dozen oscillations 
are seen in similar measurements for the latter systems, 
while we can only reliably identify 3-4 in Fig. 2). 

The simplest analysis of SdH oscillations is achieved 
by plotting positions of p xx minima and maxima as a 
function of 1/B. The Landau fan diagram thus ob- 
tained is shown in Fig. 3a. It is already clear from 
linear fits in the figure that charge carriers in our de- 
vice are characterized by a Berry phase of 37r. Fit pa- 
rameters are quantified in Fig. 3b, c, which show Bp 
and j3 as a function of V g . Linear fits of Bp(n s ) in Fig. 
3b and using dn s /dV g = 7.59 • 10 10 cm~ 2 determined 
from the measured low-field Hall constant yield LL de- 
generacy = 4 ± 0.1, in perfect agreement with the 
expected valley x spin degeneracy. Similarly, the average 
Berry phase in Fig. 3c is j3 = 2ir ■ (1.5 ± 0.1). 

We can pursue an analysis of the measured SdHO data 
further by subtracting the smooth component of magne- 
toresistance, which is mainly governed by weak localiza- 
tion, and fitting the oscillating part with the ALK for- 
mula (27l | . Varying only 3 parameters, Bp. r q , and m Cl 
we obtain high fidelity fits by imposing a physical con- 
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straint, r q < 200 is. The obtained Bp values are within 
the symbol size of those shown in Fig. 3c. The cyclotron 
mass of charge carriers obtained in such fitting is shown in 
Fig. 3d. It shows marked variation with carrier density, 
increasing towards the CNP, n s = 0, which is opposite 
to that observed in monolayer graphene 0, A fit to 

the divergent behaviour, m c — g^3(7r?i s )~ 1 / 2 , expected 
for I = 3 quasiparticles in the ABC trilayer, is shown by 
the dotted line in Fig. 3d. It yields 71 « (0.5 ± 0.1) 
cV (using v « 10 6 m/s), in very good agreement with 
the 71 value expected for few-layer graphene [26]. This 
agreement provides an important credibility check for our 
analysis, and for r q shown in the inset of Fig. 3c. In fact, 
it is quite remarkable that we are able to obtain the in- 
terlayer hopping parameter from the SdHO data. 

Further support for the existence of I = 3 charge car- 
riers is provided by the unusual QHE, which develops 
in our sample at T ps 0.35 K in fields above sa 15 T, 
Figure 4. Plateaux of Hall conductivity are observed at 
a xy — ±6e 2 /h, ±l0e 2 /h, with a step of I2e 2 /h be- 
tween the hole and electron gases across the N = LL, 
confirming its 12-fold degeneracy. This, perhaps, is most 
clearly seen in Fig 4c, which shows the magnetic field 
dependence of the longitudinal, p xx , and Hall, p xy , resis- 
tivities, for V g = 70 V (n s w 4- 10 12 cm -2 ). The minima 
of p xx occur near the center of p xy (cr xy ) plateaux, at LL 
filling factors v = n s &o/B = 6, 10, 14, as indicated by 
arrows. The behaviour near the N = LL is in stark 
contrast to that in the monolayer and bilayer graphene, 
where the first QHE plateau develops at a xy — ±2e 2 /h 
(y = 2) and a xy = ±4e 2 /h (y = 4), respectively. There 
is only a weak anomaly in our sample for n s below the 
fist QHE plateau at v = 6, which can be associated with 
a xy — ±3e 2 //i, v = 3, and is probably an indication of [14 



the developing spin-splitting, Fig. 4d. 

To conclude, we have discovered that new, I — 3 chiral 
fermions indeed exist in the ABC trilayer graphene and 
govern properties of realistic graphene devices, so they 
can be detected in experiment. These quasiparticles ac- 
cumulate a Berry phase of 37r along cyclotron trajectories 
and acquire unusual LL quantization in magnetic field, 

(vhy/2eB/(hc)) 3 

e = ±- ^ '— y n\n — l)(n — 2), n integer [4j- 

6], and therefore are revealed in magnetotransport mea- 
surements, such as we presented here. Not only our re- 
sults provide experimental validation for the large body 
of important recent theoretical work ji- 11|, but, perhaps 
more importantly, they greatly extend the perceived ex- 
perimental limits, uncovering new exciting possibilities 
for future studies. Our findings are also significant be- 



cause they establish the experimental feasibility of de- 
ploying the unusual pro pert ies of the I = 3 chiral fermions 
predicted theoretically [10], such as super-screening and 
band gap tuning, in realistic, large-area graphene devices. 
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